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|> : ABSTRACT 
O 

o 

^ ' We report the discovery of a substellar-mass companion to the KO-giant HD 

17092 with the Hobby-Eberly Telescope. In the absence of any correlation of 
the observed 360-day periodicity with the standard indicators of stellar activity, 
^ ■ the observed radial velocity variations are most plausibly explained in terms of a 

Keplerian motion of a planetary-mass body around the star. With the estimated 
■ stellar mass of 2.3M0, the minimum mass of the planet is 4.6Mj. The planet's 

I . orbit is characterized by a mild eccentricity of e=0.17 and a semi-major axis of 

^ ■ 1.3 AU. This is the tenth published detection of a planetary companion around 

^ ■ a red giant star. Such discoveries add to our understanding of planet formation 

around intermediate-mass stars and they provide dynamical information on the 
evolution of planetary systems around post-main sequence stars. 

> 

. Subject headings: planetary systems-stars: individual (HD 17092) 

ON 
p 

o 
i> 
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After a decade of searches for planets around Sun-like stars, it has become apparent that 
^ , to achieve a satisfying level of understanding of planet formation and evolution, the surveys 
^ ! have to be extended to other types of stars. So far, the most successful of these have been 
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searches for planets around low-mass dwarfs, which, among other goals, have been driven 



by anticipation that Ear th-mass planets can be found in their habitable zones (IMarcy et al. 



20051 : iMayor et al.ll2005l ). Surveys of white dwarfs, which probe ancient planetary systems - 
survivors of the evolution of their par ent stars - exemplify an extension of planet searches to 
the endpoint of stellar evolution (e.g. iKepler et at .1120051 ). Finally, searches for neutron star 
planets can provide information on planets around massive stars (IThorsett fc Dewey||l993r) 
and on planet formatio n in extreme, post-supernova environments (IKonacki fc Wolszczan 
2003l : IWang et al.lbooeh . 



Yet another, so far meagerly explored area of the extrasolar planetary research involves 
searches for planets around giant stars. More than a decade ago, precision radial velocity 
(RV) studies have established t hat GK-giant stars exhibit RV variations ranging from days 
to many hundreds of days (e.g. IWalker et al.lll989l : iHatzes &: Cochranlll993l . Il994l ). Enough 
observational evidence has been accumulated to identify three distinct sources of this vari- 
ability, namely stellar pulsations, surface activity and a presence of substellar companions. 
As Doppler searches for planets around main sequence (MS) stars become inefficient for 
spectral types earlier than F6-F8, because of paucity of spectral features and their rota- 
tional broadening, extending studies of planetary system formation and evolution to stellar 
masses substantially larger than 1 Mq is observationally difficult. A potentially very effi- 
cient, indirect way to remove this difficulty is to conduct surveys of post-MS giants. These 
evolved stars have cool atmospheres and many narrow spectral lines, which can be utilized 
in RV measurements to give an adequate precision level (<10 m s^^). Discoveries of planets 
around post-MS giants, in number s comparable to the current statistics of planets around 
MS-dwarfs (e.g. iButler et al.l 120061 ) , will most certainly provide the much needed informa- 
tion on planet formation around intermediate mass MS-progenitors (> I.5M0) and they will 
create an experime ntal basis with which to study dynamics of planetary systems orbiting 
evolving stars (e.g. iDuncan fc Lissaued Il998l ). Sufficiently large surveys of post-MS giants 
should furnish enough planet detections to meaningfully address the problem of a long-term 
survival of planetary systems around stars that are off the MS and on their way to the final 
white dwarf stage. 



I n order to addres s the above issues, we have joined the existing surveys (e.g. lHatzes et al 



20061 : ISato et al.l 120071 . and references therein) with our own long-term project to search for 
planets around evolved stars with the 9.2-m Hobby-Eberly Telescope and its High Resolution 
Spectrograph. The sample of stars we have been monitoring since early 2004 is composed of 
two groups, approxi mately equal in nuin bers. The first one falls in the "clump giant" region 
of the HR-diagram (j Jimenez et al.lll998l ). which contains stars of various masses over a range 
of evolutionary stages. The second group comprises stars, which have recently left the MS 
and are located ~1.5 mag above it. Generally, all our targets, a total of >900 GK-giants 
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brighter than ~11 mag, occupy the area in the HR-diagram, which is approximately defined 
by the MS, the instabihty strip, and the coronal dividing line. If the frequency of occurence 
of planets around MS-progenitors of GK-giants is similar to that of planets around solar-type 
stars, our survey should detect 50-100 planets and planetary systems, which, together with 
the detections from similar projects, will provide a firm basis for studies of planetary system 
formation and evolution around > IMq stars. 

In this paper, we describe our survey and the detection of a planetary-mass companion 
to the KO-giant HD 17092. Details of the observing procedure and survey strategy are given 
in Section 2, followed by a description of the basic properties of HD 17092 in Section 3. 
The analysis of radial velocity and line bisector measurements, including a discussion of the 
HIPPARCOS photometry of the star, is given in Section 4. Finally, our results are discussed 
and summarized in Section 5. 



Observations 



Observations were made between 2004, January and 2007, March, with the Hobby- 
Eberly Telesc ope (HEX] (IRamsey et al.lll998l ) equipped with the High Reso lution Spectro- 
graph (HRS) (ITu11||1998| ) in the queue scheduled mode (IShetrone et al.l 120071 ). The spectro- 
graph was used in the R=60,000 resolution mode with a gas cell {I2) inserted into the optical 
path, and it was fed with a 2 arcsec fiber. The observing scheme followed st andard practices 
imple mented in precision radial velocity measurements with the iodine cell (IMarcy fc Butler 
I992I ). The spectra consisted of 46 echelle orders recorded on the "blue" CCD chip (407.6- 
592 nm) and 24 orders on the "red" one (602-783.8 nm). The spectral data used for RV 
measurements were extracted from the 17 orders, which cover the 505 to 592 nm range of 
the I2 cell spectrum. 

Originally, HD 17092 had been observed as part of the astrometric reference star selec- 
tion pro gram related to a search for terrestrial-ma ss planets with the Space Interferometry 
Mission (iGelino et al.ll2005l : iNiedzielski et al.ll2005l ). The star was added to the list of candi- 
dates for substellar companions, when it became clear that it exhibited RV variations, which 
disqualified it as a potential astrometric reference standard. 

The observing strategy is illustrated in Fig. la. Measurements of a particular target 
star begin with 2-3 exposures, typically 3-6 months apart, to check for any RV variability 
exceeding a 30-50 m s~^ threshold. If a significant variability is detected, the star is scheduled 
for more frequent observations, and, if the RV variability is confirmed, it becomes part of 
the high priority list. 
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We have collected radial velocity measurements of HD 17092 at 59 epochs. Typically, 
the signal-to-noise ratio per resolution element in the spectra was 200-250 at 594 nm in 
3-8 minutes, depending on the atmospheric conditions. The basic data reduction was per- 
formed using standard IRAF0 scripts. Rad ial velocities were measured by means of the 
commonly used I2 cell calibration technique (IButler et al.lll996l ). A template spectrum was 
constructed from a high-resolution Fourier Transform Spectrometer (FTS) I2 spectrum and 
a high signal-to-noise stellar spectrum measured without the I2 cell. Doppler shifts were 
derived from least-square fits of template spectra to stellar spectra with the imprinted I2 
absorption lines. The resultant radial velocity measurement for each epoch was derived as a 
mean value of the independent determinations from the 17 usable echelle orders. The corre- 
sponding uncertainties of these measurements were calculated assuming that errors obeyed 
the Student's t-distribution and they typically fell in a 4-5 m s~ ^ range a t Icr-le vel. Radial 
velocities were referred to the Solar System barycenter using the IStumpffl (jl980l ) algorithm. 



The star 



HD 17092 (BD+49 767) was classified as a KO-star by ICannon fc Pickerind fll924h . 
These authors also measured its photographic and pho tovisual magnitud es to be Ptg=8'".8 
and Ptm=7'".8, respectively. The Tycho-2 catalogue (iHogg et al.l l2000l ) lists the values of 
Br=9'^.374±0.019 and V^=7'".875 ±0.011 der ived from the HIPPARCOS observations of 
the star. In the Tycho-1 catalogue JESAlll997h .values of V=7™.73, B-V=l"".247±0.014 are 
also given. In the Tycho star mapper experiment, a trigonometric parallax of HD 17092 was 
measured as 9.2±5.5 mas. 

A detailed inspection of our spectra in the range of 515-520 nm (comparison of intensities 
and widths of Mgl b triplet lines at 516.7, 517.2 and 518.4 nm) suggests that the star is a 
giant . The absolute visual magnitude of HD 17092 is My=l'".76, assuming (B-V)o=l.'"00 
after ISchmidt-Kalerl (119821 ). the above Tycho parallax, and R=3.1. This makes the sta r 
almost one magnitude fainter compared to a typical Mv'(K0III)=0'".8 (jSchmidt-Kalerlll982l ) . 
In what follows, we will assume that HD 17092 is a typical giant with My=0'".8, as this 
value remains consistent with the Tycho parallax uncertainty. 



The atmospheric parameters of HD 17092 were recently estimated by lZielinski fc Niedzielski 



(120071 ) through the analysis of 230 Fel and 11 Fell lines in the optical spectra. The results. 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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Te//=4650K, log;(K )=3.0, and [Fe/H]=0 . 18, ar e very close to rough estimates obtained from 
the calibrations bvlStraizvs fc Kuriliend (Il98ll ). which give logTe//=3.681, and log(g)=2.89. 



have also estimated the mass of HD 17092 as M /M,t,=2.3±0.3 



by com paring the star's position in the HR diagram with evolutionary tracks of lGirardi et al. 



( I2OOOI ) given the above metalicity. In absence of a direct measurement of the radius of 
HD 17092, we assume that it is similar t o the radii of stars with [Fe/H]> -0.5 and adopt 
R/R0=1O.1±4.2 after IaIouso et all J2000h . 



The projected rotational velocity of HP 1709 2, v sini <1 km s^^, was estimated using 
the cross-correlation method (IBenz &: Mayorlll984j ). From this value and the adopted stellar 
radius we have obtained an estimate of the ro tation period of HP 170 92, P^ot ~505 days. 
This value appears to be typical for KO giants (jde Medeiros et aL 1996 ). 



4. Data analysis 
4.1. Modeling of the companion orbit 



Radial velocity variations of HP 17092 over a 3-year period are shown in Fig. la, to- 
gether with the best-fit model of a Keplerian orbit. The best-fit parameters of the orbit 
and their Monte Carlo estimated uncertainties are listed in Table 1. The residuals shown 
in Fig. lb are characterized by the rms value of ~16 m . If the observed RV vari- 
ations are indeed caused by an orbiting companion, it moves in a ~360-day, moderately 
eccentric orbit, with a semi-major axis of 1.3 A.U., and has a miniumum mass of m.2 sin 
?=4.6Mj, for the assumed stellar mass of 2.3 Mq. This mass indicates a planetary origin 
of the object over a range of possible values of sin i extending beyond the median incli- 
nation of i = 60° for randomly oriented orbits. The reduced x^=10 for the fit suggests a 
presence of additional variati ons in the data, w hich have also been reported for other planet 
detections around gia r its fe.g. iFrink et al.ll2002r) and are not uncommon for this type of stars 
dSetiawan et~aD 120041 : iHatzes et allboosh . Our work also confirms that red giants exhibit a 
RV scatt er at an average level of ~20 m s~^, as the result of a stochastic intrinsic stellar 
activity (INiedzielski fc WolszczanI 120071 ). To account for these variations, we have adopted 



a conservative error estimation procedure by adding in quadrature a constant 15 m s~^ er- 
ror to the RV measurement uncertainties, before performing a least-squares fit of the orbit. 
This approach resulted in parameter error estimates (Table 2), which realistically absorb any 
leftover, unmodeled RV variations in the data. Their nature will be further discussed, when 
more data become available. 
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4.2. Line bisector and curvature analysis 



Precision radial velocity measurements may be significantly affected by phenomena that 
are not related to stellar reflex motion caused by the presence of an orbiting planet. Changes 
in line shapes arising from motions in the stellar atmosphere, related to non-radial pulsations 
or inhomogeneous convection and/or spots combined with rotation, or distortions induced 
by light contamination by an unseen stellar companion, can mimic low-level radial velocity 
variations. Therefore, especially in the case of giant stars, it is important to verify whether 
the observed radial velocity variations are real, or if they are possibly generated by changes 
in the symmetry of spectral lines due to the above effects. 

The basic tool to study the origin of RV variations derived from stellar spectra is the 
analysis of shapes of the spectral lines using the line bisector technique (lGraylll983l . l2005l ) . It 
has been extensivel y used in the process of con firmation of a planetary origin of RV variations 

1998a|]b[) and it has become a mandatory part of the 



in solar-type stars (IHatzes et al 



19971. 



2005 



analysis of RV data from giants (ISetiawan et al.l l2003l : ISetiawan et al.l l2005l : iHatzes et al. 



Reffert et al.ll2006l : ISato et al.l 120071 ). To examine line profile variations in HD 17092 



we have selected several lines of a moderate intensity, which were free of blends and were 
located close to the centers of echelle orders. Because the HRS spectra extend far beyond 
the range occupied by the I2 lines, we were able to measure line profiles in the same spectra 
which were used for radial velocity determination. 

In order to avoid any contamination by the I2 spectrum, we have selected lines in 
the wavelength range redwards of 660 nm, which also meant that we could not use the 
standard bisector line of Fe I at 625.256 nm. Instead, we have used three lines, Ni I 664.638, 
Ni I 676.784, and Cr I 6 63.003, whi c h wer e analyzed in detail to search for any possible 
systematic effects. In fact. iDall et al.l (120061 ) has found that the Ni I 664.3683 nm line shows 
larger bisector variations than the Fe I 625.256 nm line. Also, the other two lines adopted 
for this analysis show well defined bisectors. We have measured two line parameters, the line 
bisector span and its curvature, under the assumption that the mean bisector and velocity 
span of the three lines reflect the upper limit of a possible atmospheric activity of the star. 
Uncertainties in the derived values of the bisector span and curvature were estimated as 
standard deviations from the mean. As shown in Fig. 2, both parameters, when plotted as 
a function of orbital phase, are not correlated with radial velocities. 
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4.3. Photometric variabilty 

The HIPPARCOS star mapper (TYCHO) has made 124 photometric observations of HD 
17092 in Vt and Bt between JD 2447915.85944 and 2449039.86213, about 12 years before 
the beginning; o f our survey. The measured scatter in Vt (defined as 0.5(Vt(85)-Vt(15)), 



see 



ESAl ( 119971 ) for details) was 0.092 mag without any sign of systematic variability. The 
observed scatter is consistent with the precision of Tycho photometry for a given stellar 
magnitude range. 

We have performed a deeper search for any possible periodicities in the Tycho photom- 
etry of the star by computing a Lomb-Scargle periodogram of these data. No significant 
peaks were found in the spectrum above the false positive probability of 0.5 (Fig. 3). In 
particular, no excess fluctuation power is present at and around the 360-day period detected 
in the RV data. 



Discussion and conclusions 



In this paper, we have presented a compelling evidence that the KO-giant, HD 17092, 
exhibits a strictly periodic radial velocity variation with the period of 360±2 days over a 1200- 
day span of observations. When interpreted in terms of a Keplerian motion, this periodicity 
indicates the presence of a sub-stellar companion with a minimum mass of 4.6 Mj, in a 
moderately eccentric orbit (e=0.17±0.05), 1.3 A.U. away of the star. As the long period 
RV variations in red giants may also be related t o a combination of effects including stellar 
rotation, activity, and non-radial pulsations (e.g. iHatzes et al.ll2006l ). we have analyzed the 
photome tric data and the b ehavior of line bisectors of the star, following the established 
practice (iQueloz et al.ll200ll ). 



The magnitude of a photometric variability is related to a fraction of the stellar surface 
covered by spots. The Tycho photometry, although not very precise, gives an upper limit 
to this variation of 0.011 magnitud e in Vr which translates to ~1% of the stellar surface 
covered by spots. As discussed by iHatzed (120021 ). with this spot coverage and the star's 
rotational velocity of ~1 km s~^, it is impossible to produce RV variations of the observed 
amplitude. In fact, their maximum amplitude would be on the order of a few m s"^, which 
is comparable to the precision of our RV measurements. 



Hatzed (120021 ) has published estimates of the bisector velocity span (BVS) as a function 



of vsini and fraction of spots. For HD 17092, we derive from their formula the BVS of a few 
m s~^, which is consistent with our own BVS determination. A detailed analysis of both line 
bisectors and bisector velocity span shows that there is no relationship between line profile 
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variability and the observed RV changes. In addition, hne profile variations, measured with 
the same spectra as the ones used to obtain RVs, show no correlation with RV. 

A lack of both line profile variations and photometric variability eliminates nonradial 
pulsations and surface activity related to stellar rotation as a possible cause of the observed 
RV variations in HD 17092. Another argument against a rotation-forced RV variation is 
the repeatability of the observed RV changes and their strictly periodic character. Spots, as 
we understand them, appear and disappear on the stellar surface and the number of spots, 
as well as their location, vary from one cycle to another. Therefore, a resulting variability 
is unlikely to maintain the same pattern over many cycles. In the case of HD 17092 the 
observed 3 cycles of the 360-day period are consistently fitted with a single Keplerian orbit. 
Consequently, the most plausible explanation of our data is the presence of a sub-stellar 
mass companion around the star. 

A sufficiently large number of planet detections by high precision RV surveys of GK- 
giants will make them efficient tools with which to study planet formation around intermediate- 
mass stars and the dynamical evolution of planets induced by a post-MS evolution of their 
parent stars. Current constraints on a stellar mass dependence of the disk mass and the 
timescale of depletion of its gas and dust components come fro m studies of disks around 
young stars. For example , in addition to the pr evious work (e.g. iHaisch et al.ll200ll ). recent 



Spitzer observations (e.g. ICarpenter et al.ll2006l ) appear to confirm that disks around inter- 



mediate and high-mass stars have lifetimes significantly shorter than 5 Myr. These results 
have direct consequences for the competing theo r ies of giant planet formation, because the 
core accumulation scenarios (jPoUack et al.l Il996l : IWuchterl et al.l l2000l : lAlibert et al.l |200J) 
require at l east a few rnillion years for a c ore to form, whereas planet formation from a disk 
instability (IBossI 119971 : iMayer et al.ll2002l ) can be very short. Clearly, the searches for plan- 
ets around giant stars have a unique capability to provide the statistics which are needed 
to decisively constrain the efficiency of planet formation as a function of stellar mass and 
chemical composition. 

It is quite reasonable to expect that the giant star surveys will have the potential to 
veri fy predictions of the post-MS orbital evolution which emerge from numeric al simulations 
(e.g. iRasio &: Fordlll996l : iDuncan &: Lissauerlll998l : iDebes &: Sigurdssonll2002l ). In principle. 



each detection of a planet around a giant represents a snapshot of the dynamical evolution 
of orbits around a particular, evolving star. Given a sufficient number of planet detections it 
should be possible to constrain the principal evolutionary scenarios and obtain a consistent 
picture of dynamical changes in planetary systems in response to the evolution of their parent 
stars. An excellent example of the a nalysis of planetar y orbits, which would be applicable 
in this context has been presented by iFord et al.l (120051 ). 
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HD 17092b is the tenth pubhshed discovery of a planet around a giant star. It is one of 
the most distant and metal-rich giants known to host a planetary system. Further discoveries 
of planets around giant stars will help to develop sufficient statistics to meaningfully address 
the important questions of planet formation around intermediate-mass stars and the long- 
term evolution and survival of planetary systems. 
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Fig. 1. — (a) Radial velocities (filled circles) and the best-fit orbit (solid line) for HD 17092. 
(b) Residuals from the best-fit of a Keplerian orbit to data. The vertical bar indicates an 
additional error added in quadrature to the formal uncertainties of radial velocity measure- 
ments. 



-14- 




Fig. 2. — The mean bisector velocity span (a) and curvature (b) of HD 17092, binned as a 
function of orbital phase at the 0.1 interval, and compared to the measured radial velocities 
of the star (c) . 
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Fig. 3. — A periodogram of the HIPPARCOS photometric measurements of HD 17092. The 
vertical arrow marks the frequency corresponding to the 360-day period of the observed 
radial velocity variations. 



Table 1. Stellar parameters of HD 17092 



Parameter Value 



V 

B-V 

Spectral type 

log g 

[Fe/H] 

M* 



7.73 

1.247±0.014 
KOIII 
4650±35 
3.0±0.12 
0.18 it 0.08 
2.3 ±O.3M0 
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Table 2. Measured and derived orbital parameters of HD 17092b 



Parameter 


Value 


P (days) 


359.9 ± 2.4 


To (MJD) 


52969.5 ± 12.3 


K (m s-l) 


82.4 ± 3.2 


aisin i (AU) 


0.0027 ± 0.0002 


e 


0.166 ± 0.052 


u) (dcg) 


347.4 ± 13.4 


f(m) (Mq) 


2.001x10-8 ± 1.1x10-9 


m2sin i (Mj) 


4.6 ± 0.3 


a.2 (AU) 


1.29 ± 0.05 



